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I. INTRODUCTION 

Using the AdS/CFT correspondence [l[ it has been shown that if an abelian symmetry has been broken outside of 
a black hole in AdS space, scalar hair can form creating a holographic superconductor in the dual CFT0-11]. In the 
last few years it has been seen that these systems exhibit several characteristics seen in real world strongly coupled 
superconductors and seem very promising. The AdS/CFT correspondence has also been applied to other areas of 
condensed matter physics 043 > some reviews are [Iol - [T2j . 

Studying the conductivity of holographic superconductors, Horowitz and Roberts [1, [H| noted that at the 
Breitenlohner-Freedman (BF) bound of the scalar hair of the dual black hole, quasinormal modes appeared to 
move toward the real axis at low temperatures. Once the back reaction to the metric was included, they showed that 
the quasinormal modes never became normal. Their number was determined by the height of an effective potential 
associated with the wave equation of an electromagnetic perturbation. The study was done numerically which limited 
the ability to go to very low temperatures. 

Using the analytic tools developed recently in [l5j . we explore the low temperature regime of holographic super- 
conductors near the BF bound. We find that in the probe limit, when the BF bound is saturated, the condensate 
diverges at low temperatures as | In where S is a constant that we compute and depends on the dimension of 
spacetime. This is a very mild divergence which explains why it was not detected in earlier numerical studies |4[. It 
signals the breakdown of the probe limit at very low temperatures, and by including corrections to the probe limit 
we determine when this occurs. When back reaction to the metric is included, the effective potential associated with 
the electromagnetic perturbation that determines the conductivity has a finite height. This results in a finite number 
of quasinormal modes. As one approaches the probe limit, the height of the potential increases with an attendant 
increase in the number of modes. The latter approach the real axis as the temperature is lowered at frequencies 
that we compute analytically. In 2+1 (3+1) dimensions they are given in terms of the zeroes of the Airy (Gamma) 
function. 

We emphasize that, even though we probe the low temperature regime, we do not have access to the zero temperature 
state. This is because for a given charge of the scalar field, there is a lower bound on the temperature (for scaling 
dimensions A < f ) below which the probe approximation breaks down. Thus, even though it is often possible to take 
a "zero temperature" limit of our expressions, the limit itself is unphysical. It is, however, useful for computational 
purposes, as physical low-temperature systems are close to it, as we shall show. 

Our paper is organized as follows. In section ITT1 we review the field equations. In section Hill we calculate the critical 
temperature. In section HVl we discuss the low temperature regime in the probe limit. We also discuss the validity 
of the probe limit, and find where it breaks down for a given value of the charge q. In section [V] we calculate the 
conductivity at the BF bound at low temperatures. Finally, section fVTl contains our concluding remarks. 



II. FIELD EQUATIONS 



We are interested in the dynamics of a scalar field of mass m and electric charge q coupled to a U(l) vector potential 
in the backgound of a d + 1— dimensional AdS black hole. The action is 



S = / d d+1 J= 



<) 



R + d{d-\)/l 2 1 



IQnG 



(1) 



where F = dA. We shall adopt units in which 1 = 1. 

To find a solution of the field equations, consider the metric ansatz 



ds 2 = \ 



-f(z)e~ x ^dt 2 + dx 2 + 



dz 1 

W) 



(2) 



where x S representing an AdS black hole of planar horizon. The AdS boundary is at z = 0. We shall choose 

units so that the horizon is at z = 1, therefore we require /(l) = 0. This is possible because of scaling symmetries 
of the system and can be done without loss of generality as long as one is careful to only consider physical quantities 
which are scale invariant. 
The Hawking temperature is 



T 



An 



-x(i)/2 



(3) 



3 



Assuming that the scalar held is a real function ^(z) and the potential is an electrostatic scalar potential, A — $(z)dt, 
the field equations are Q 



/' x' 



m 



z 2 f 



* = 



x_ _ 
2 

d- 1 , 



d- 3 



2o 2 * 2 
z 2 / 
zq 2 ^ 2 ^ 2 



L^ 2 

2 



1/' d(d-l)/-l m 2 * 2 q 2 * 2 $ 2 e x 



2z 2 



2/ 



(4) 



where prime denotes differentiation with respect to z, to be solved in the interval (0, 1), where z = 1 is the horizon 
and z = is the boundary. 

We are interested in solving the system of non-linear equations Q in the limit of large q (probe limit). To this end, 
we shall expand the fields as series in 1/q as follows: 



1 



q- 

1 



/ = /o + /i- + - 

l 

X = Xo + Xi^r + 



(5) 



and consider the zeroth order system (q —> oo) first and then discuss the addition of first-order (0(1/ q 2 )) corrections 
in order to obtain a physically sensible system. 

Near the boundary (z — > 0), we have / — > 1, % — > and so approximately 



$ fa /i — pz 



d-2 



where 



2 V 4 



(6) 



(7) 



While a linear combination of asymptotics is allowed by the field equations, it turns out that any such combination is 
unstable [20| . However, if the horizon has negative curvature, such linear combinations lead to stable configurations 
in certain cases [2lj . 

Thus, the system is labeled uniquely by the dimension A = A±. The mass of the scalar field is bounded from below 
by the BF bound [HI], m 2 > — 4-, and there appears to be a quantum phase transition at m 2 — 0. There is also a 
unitarity bound that requires A > 
Demanding at the horizon 



d-2 



$(1) = 0, (8) 

(gauge choice ensuring that A = Qdt is regular at the horizon Q), [i is interpreted as the chemical potential of the 
dual theory on the boundary, p is the charge density on the boundary and the leading coefficient in the expansion of 
the scalar yields vacuum expectation values of operators of dimension A± , 



V2*<= 



(9) 



The field equations admit non- vanishing solutions for the scalar below a critical temperature T c where these operators 
condense. In view of ([6]) and (J9j), it is convenient to define 



*(z) 



f 



V2q 



b /X z^F(z) , b=(qO A 



(10) 



with F(0) = 1. 
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III. THE CRITICAL TEMPERATURE 



Above the critical temperature, — and the field equations are solved by the AdS Reissner-Nordstrom black hole 
with flat horizon, 



/(*) = 1 - (l + z d + V-3f z ^-V , x{z) = o , = p (1 - z d-*) 

whose Hawking temperature is 



(11) 



T 



4tt 



(rf-2)V 
4d 



(12) 



Right at the critical temperature, \E» obeys the scalar field equation (jj) in the Reissner-Nordstrom background (fTTj) 
with p — p c . Thus -F (eq. (jTUJ) ) at T = T c obeys the field equation 



F" 



f 2A + l-rf 

— H 

/ 



■ A (d-A)(l-/) + z/' 2 2 (l-z d - 2 ) 2 

A — + q p c 



z 2 f 



p 



F = 



(13) 



For a given g, /9 C is an eigenvalue which is determined by solving this equation for F subject to the boundary condition 
at the AdS boundary F(0) = 1. We also demand that at the horizon F(l) be finite, and that there be no contribution 



of the other solution (behaving as F(z) 



asz-> 0). The latter condition implies that F has a Taylor expansion 



around z — with the properties F(0) = 1, F'(0) = (as can be easily deduced from ((T 
To solve eq. (fT3| for large q (probe limit), use the expansion §5§ to write 



F — Fq + Fi - 



1 



1 



P 



1 

Po + Pi -« + 



<r q 

Then at zeroth order (q — > oo limit), the background (TTTT) turns into an AdS Schwarzschild black hole, so 

fo(z) = l-z d , Xo = 
and we obtain the field equation at the critical temperature 



F'' 



1 



- 2A - 1 



„d-2 



2 (l_ 2 d-2)2 
F ° + ~ ^ ° = P ° C (1-Z d ) 2 ° 



1 - Z d 

The eigenvalue Pq c minimizes the expression 

2 _ Jo 1 dzz 2A - d+1 {(l - z d )[F^z)] 2 + A 2 z rf - 2 [F (z)] 2 } 

We can estimate the eigenvalue by substituting the trial function 

F =F a (z) = 1 -az^ 1 
which obeys the boundary conditions F a (0) = 1, i^(0) = and F Q (1) is finite. 



(14) 



(15) 



(16) 



(17) 



(18) 



For A — | and rf = 3, 4, we obtain, respectively, 



p 2 0c « 6.3, 4.2 , T c a 0.15^, 0.2(gp) 



1/3 



(19) 



in very good agreement with the exact T c = 0.15^/gp, 0.25(qp) 1 ^ 3 . In fig. [T] we extend the comparison to the entire 
range of the scaling dimension A for d = 3, 4 demonstrating the accuracy of the estimate (1171) with the trial function 
(JTSJ) for the critical temperature. 
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FIG. 1. The critical temperature T c vs the scaling dimension A for d = 3 (left panel) and d — 4 (right panel). Data points 
represent exact values; solid line is obtained by minimizing (|f 7p with the trial function JT8J. 



IV. LOW TEMPERATURES 



Next we consider the low temperature regime. Because we are working in units in which the radius of the horizon 
is fixed (z — 1), even at low temperatures, T is bounded. Therefore, p becomes large. Also the condensate diverges in 
these units for the same reason, and so does b (eq. (|10[) ). We are interested in calculating scale-invariant quantities, 
such as the energy gap 



l/A 



!/(<*-!) 



(20) 



We shall consider the probe limit (q — > oo) and then discuss how first-order corrections (in a 1/g expansion) can be 
added to obtain a physically meaningful system. 

In the probe limit, as we lower the temperature, Fq(z) (eq. (|14[) ) does not have a smooth limit as b —¥ oo for 



A = A+ > §. In this case, care should be exercised in taking this limit which corresponds to the zero temperature 
limit (ground state) of the system. 

We shall concentrate on the case A = A_ < i in which Fq{z) has a smooth limit as b — > oo, however the limit 
system does not correspond to the zero temperature limit and is not physical. This is because the probe approximation 
breaks down as b — > oo. As we shall show, the probe expansion is valid for b < q x / A (eq. Q5ip). So for a given charge 
q ^> 1, we define the low temperature regime by 



1 Z,b<q 



l/A 



(21) 



We need to solve the system of zeroth-order equations, 



l-z° 



1 - 2A 



. A 2 z d - 2 



(I - z d ) 



z 1 — z a 



(22) 



Let us first discuss the b — > oo limit. Even though it is unphysical, it is well-defined and will be useful for our analytic 
calculations. 

It is evident from the field equation for $o (and can be easily confirmed numerically for arbitrary regular functions 
F (z)) that $o->0 for z>l/b (with b > 1). Then for z > 1/6, we obtain 



(23) 



6 



which is regular at the horizon. As b —> oo, this becomes valid in the entire interval, because 1/6 — > 0. We deduce 
the b — > oo function 



r»(i-4 



r(i 



A A 



P • 1 • 1 - z '- 

M d' d' ' 



(24) 



where we used Fo(0) = 1 and standard hypergeometric identities. However, this does not represent the ground state of 
a physical system, because as b — > oo, the condition for a physical system ([21]) is violated. Nevertheless, the expression 
(|2"41 is useful for computational purposes. 

Next, we wish to solve the equations in the low temperature regime (|2"Tj) . To this end, we shall use iteration as 
follows, 



F, 



(n+l) 







" 1 
+ - 

z 



1 - z d 



1 - 2A 



F, 



(n+l)' i A z jp(n+l) 



l-z~ F ° 



ii rA(»+l)l2 F (n) 



(n+l)// d-^3^(„+l)/ 52A-,2(A 1) rF ( n)l2 ^(„ +1 ) 



1 - Z fl 



[f}"fC = o 



starting with 



We defined 



F (0) (z) = l , <^ 0) (z) = 



$oW = AW , *o(0) = l 



(25) 



(26) 



(27) 



where /x is the chemical potential. 

At the nth step, we obtain the electrostatic potential $q™ +1 ' ) from the second equation in (|25|) and subsequently the 
scalar field F^ n+1 ^ from the first equation in (1231) . To solve the latter, notice that there are three boundary conditions, 
two at z = (i*o (0) = 1 and absence of the unwanted alternate behavior Fq ~ z d ~ 2A ), and one at the horizon, 
z = 1 (finiteness of F (n+1) ). They are needed to determine the solution as well as the eigenvalue \i. Using the two 
boundary conditions at z = 0, we obtain the solution in terms of the chemical potential, 



F, 



(n+l) 







(z)=T 1 (z) 



dz' 



/\2A + l-d 



j^0[^ n+ V)] 2 ^o (n V) 



where 



Fi(z) = F 



A A_2A 

d' d ' d 



o l-(^') c 
, ^(/) 2A+1 - d J- 1 (.')[*^ +1) (^)] 2 ^o (n) (^) 



o l-(*0 



V t2-2A 



-F 2 (z) 



d- 2A 



F l--,l--;2- 



2A 



(28) 



(29) 



The third boundary condition, regularity at z = 1, then fixes the chemical potential /i. 
For n = 0, we obtain the electrostatic potential 



r(i/)(2A)= 



:(&*) 



A 



^ "l A 



d- 2 
2A 



(30) 



where we imposed the boundary condition (jHJ. Notice that the second Bessel function has an exponentially small 
coefficient, 0(~ e" 2b / A ), and can be neglected at low temperatures. 
The charge density is found by using © and (fl"4"]) to be 



fe d-2 



(2A) 2 ^ 



For the scalar field we obtain eq. 



with n = 0, 



F (1) (z) = J-r(z) 



1 + M 



^7w(^) 2A+1 - d ^(,')[$^(/)] 2 



l-(z') 



f z dz' 



(31) 

'^ +1 - d Fi{z')[^\z')f 
(32) 
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FIG. 2. The field F (eq. CGI) for A = 1.2 (left panel), 1.4 (middle panel), 1.5 (right panel) and d = 3. Solid curves are 
first-order analytic expression (|32[1 . and dashed curves are exact numerical results (almost indistinguishable) at T/T c ~ 0.1. 



satisfying the correct boundary conditions at z = 0. At the horizon (z — 1), we have a logarithmic singularity which 
is found using 



Near the horizon, we deduce 
where 



r 2 (£) 



ln(l-z) , F 2 {z) 



T(2-2A) 



(d-2A)I*(l-£) 



Ax ln(l-z) 



F(2A) 



r(2 



2A ' 



dz 



„2A + l-d 



(d-2A)r 2 (i- 



ln(l - z) 



1 



Demanding regularity at the horizon, we need 

'2Ai 



r 2 (f) ( 



r(2-g#) 
(d-2A)i*(i-f 



i = 1,2 



which fixes the chemical potential /i, 



r(2-^)r 2 (#) 

(d - 2A)T(2A)r 2 (l - — ) 



ai — a 2 



as advertised. 
Explicitly, 



«i 



1 (d-2)r(l-i 
5 2A+2-d ( 2 A) 2l T(i/) 



a-2 



i 0FAi-ir(^)r(^)r(^) 
6 2 (d-2A)r 2 (^)2 2 ^r(^) 



Evidently, for A < =, a 2 <C a\ for 6 <C 1, therefore 



(33) 



(34) 



(35) 



(36) 



(37) 



(38) 



fi 2 « C& : 



2A+2-d 



(d - 2A)(2A) 2t/ i»r(^)r 2 (i - § 
(d-2)r(i-^)r(2-2A)r 2 (A) 



(39) 



It is easily seen (using standard hypergeometric identities) that the low temperature expression (|32|) reduces to (|24| as 
— > oo in the entire interval [0, 1]. However, as we have already pointed out, the latter is unphysical. By taking this 
limit, we do not gain access to a zero temperature system (gound state). Nevertheless, it is useful computationally. 

Before we consider the next iterative order, we note that at finite temperature, the first-order expression (|3"2")l 
is in excellent agreement with numerical results even at T/T c ~ 0.1, which is the lowest temperature at which 
a numerical solution is available. This is shown in figs. [5] and [3] in which the corresponding curves are almost 
indistinguishable, implying that the next iterative order introduces negligible corrections to the first-order expression 
(1521 for temperatures T/T c < 0.1. 
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FIG. 4. The second order correction to the scalar field F for d = 3, A = 1.4 at b = 20 (T/T c ~ 0.1) (solid line) and b = 200 
{T/T c ~ 0.01) (dashed line). 



(2) " (2) 

We can repeat the above steps for the next iterative order to calculate F and $q . The resulting functions are 
very close to the their first-order counterparts, showing that the iteration converges rather rapidly. In fact, the second 
order quantities are subleading in 1/6. This is the case for all values of the scaling dimension A < |. In fig. HI we 

show the difference between first and second iterative order for d = 3, A = 1.4 (all other values of A are similar). 

F W 

The error, 1 tj-t, is less than 5% in the entire interval [0,1]. As the temperature decreases from T/T c ~ 0.1 to 

F o 

T/T c ~ 0.01, the error decreases to less than 0.01. To demonstrate that the error is subleading in 1/6, in fig. [5] we 
plot it at the mid-point (z = -|) as well as the horizon (z = 1) (at z = the error vanishes by design). Evidently, it 
goes to zero as 1 /6, showing that the second iterative order introduces subleading corrections at low temperature. 
For the charge density we deduce from (j3"T1) 



Po ~ b^- 1 (40) 



Using 



we finally obtain for the energy gap 



bPo'- 1 , Y~ Po (41) 



d/2-A 



E B = " Z =7[ ¥ ) (42) 



showing that the condensate diverges at low tempartures (to be precise, in the regime (|2ip). The exponent depends 
on the dimensions of the operator and spacetime. The expression for the exponent in (|42l) corrects an earlier analytic 
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FIG. 5. The second order correction to the scalar field F for d = 3, A = 1.4 as a function of temperature at the mid-point, z = 4, 
(left panel) and the horizon, 2 = 1, (right panel). The horizontal axis corresponds to the temperature range 0.01 < £ < 0.1 
with T decreasing to the right. 




FIG. 6. The parameter 7 in the low temperature expression (|42|) for the condensate vs A. Curve on left (right) is for d — 3 
(d = 4). 



result 15;j. The constant of proportionality 7 can be found analytically. It is plotted in fig. \6\vs A. As A approaches 
the BF bound, 7 — > 0, showing that the power law behavior changes, as we discuss next. 

Indeed, at the end point (BF bound), A = |, we need to exercise care. Letting A = | — e, we obtain from ([57)1 
and HI, 



and taking the limit e 



1 _ (d- 2)r(§) 

<po-5)r(i-%) 

0, we deduce at the BF bound 

^ _ d 2 ^ 
b 2 ~ 



1 



1 



(43) 



(d-2)T(j,)[lnb + (3 d + o(W)} 



(44) 



where /3d is a constant that depends on the dimension and is easily computed (e.g., for d = 3, (3$ ~ 1.75). 

Higher-order corrections are computed as before by considering the next iterative order. This introduces corrections 
that are subleading as b — » 00. Already at T/T c ^0.1, our first-order analytic results are almost indistinguishable from 
numerical results (see right panels of figs. [2] and [3]), the discrepancy being ~ 1%. Below that temperature, numerical 
results are not available. Nevertheless the error can be estimated by calculating the correction at the second iterative 
order, as before. One finds that the error tends to zero as b — ¥ 00 in the entire interval [0, 1]. 

For the charge density, we have 



Po 



b d '\lnb)-^ 2 



(45) 
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therefore the energy gap behaves as 



E g = 



{qOi 



a/A 



(In 6) ^ 



In 



(46) 



showing that the condensate diverges at the BF bound, albeit very mildly. This mild divergence was missed in earlier 
numerical studies [J]. 

The BF bound can also be approached from above. However, the calculation becomes considerably more invloved, 
because for A > d/2, as T — > 0, we have F w 1 near the boundary (z — 0), but asymptotically (z > 1/b), 
F ~ z d ~ 2A , which does not have a smooth limit as T — > 0. Therefore, we cannot apply perturbation theory and a 
different approach is called for [l5j]. For example, one can approximate F by 



F {z) 



1 



-2A 



z < a 
z > a 



(47) 



and find a by a variational method. We shall not dwell on this further here. 

Having understood the probe limit, we now turn to the first-order corections in a 1/q 2 expansion. For A < d/2, 
it is necessary to include these corrections in order to obtain a physical system at low temperatures, because in the 
q — > oo limit the condensate diverges as T — > (eqs. (l4"2l and (HHJ)- 

At first order, we obtain for the functions determining the metric, 



zf'i - dfi 



(bz 



|2A 



4(d-l) 
(bz) 2A 



A 2 /o 



2A 



(bz) 



F 2 + 2Azf F F{ ) + z 2 f Ff 



(bz) 



2A*0 



/ 2 



1 



A 



/o 2 



Fq + 2AzF F^ + z 2 F^ 



They can be solved at low temperature using our zeroth-order results above. We obtain 

A ,. 



h(z) 



-(bz) 2A [2-z d - z d ~ 2A ] + 



A(d-iy 

For the temperature, we deduce the first-order expression 



Xi (z) 



2(d-l) 



(bz 



,2A 



T = 



47T 



U2A 



1 



2d(d - 1) q 2 



(48) 



(49) 



(50) 



showing that the temperature receives a positive correction away from the probe limit. Moreover, it is now clear when 
the probe limit fails. Indeed, for the expansion in 1/q 2 to be valid, we ought to have 

b<q 1/A (51) 

justifying our earlier definition of the low energy regime f|2 1[) . 

For a given q, (l5Tj) places a lower bound on the temperature. While zero temperature is unattainable for finite q, 
the temperature can be made arbitrarily low by choosing a sufficiently large q. It follows that, even though the probe 
limit (q — ¥ oo) is not a physical system at zero temperature, its properties are a good approximation to corresponding 
properties of physical systems (of finite charge q and (low) temperature). The approximation becomes better with 
increasing q as the 1/q 2 corrections become smaller. 



CONDUCTIVITY 



Next, we calculate the low temperature conductivity at the BF bound. For explicit analytic results, we concentrate 
on two cases, d = 3 and d = 4. We shall obtain the conductivity a as a function of the rescaled frequency 



a/A 



(52) 



The function a(u>) has a well-defined limit as q — > oo (probe limit) down to low temperatures (in the regime (|21[1 ) at 
which the condensate (Oa) is large. The conductivity of physical systems can be obtained as a 1/q 2 expansion with 
the conductivity in the probe limit serving as the zeroth order term in the expansion. 
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A. Three dimensions 



The conductivity on the AdS boundary is found by applying a sinusoidal electromagnetic perturbation in the bulk 
of frequency ui obeying the wave equation 



d 2 A 
drl 



+ VA = cu 2 A 



V = 



(53) 



where A is any component of the perturbing electromagnetic potential along the boundary. Eq. (153")) is to be solved 
subject to ingoing boundary condition at the horizon 



A ~ e~ ltjr * - (1 - zY^I* 
as z — > 1 (r* — > — oo), where r* is the tortoise coordinate 



dz 

W) 



(1 — z) 3 
ln^ i- - 2V3tan 

1 — z A 



2 + z 



(54) 



(55) 



with the integration constant chosen so that the boundary is at i\ — 0. We with to solve this equation at low 
temperatures. 

Using (fTCTj) with d — 3, the wave equation reads 



d_ 

dz 



(1-z 3 ) 



dA 

dz 



b 2A z 2A - 2 F 2 (z) 



1-z 3 



A = 



To account for the boundary condition at the horizon, set 

A = (1- z)- iuJ / 3 t 



z/3 



A(z) 



(56) 



(57) 



where we included a factor e %U)Z I 3 for convenience, so that only A(z) will contribute to the conductivity. The wave 
equation becomes 



-3(1 - z 3 )A" + z [9z - 2(1 + z + z 2 )iuj] A! 



3b 2A z 2A - 2 F 2 {z) - (1 + 2z + Zz 2 )iuj - 



(3 + 2z + z 2 )(3 + z + z 2 + z 3 ) 



3(1 + z + z 2 ) 

Regularity of the wavefunction A at the horizon (z = 1) implies the boundary condition 



.4 = 



(3 - 2tw) A'(l) + b 2A F 2 (l) - 2ioj 



A(1) = 



(58) 



(59) 



At low temperatures, b 3> 1, it is convenient to rescale z — > zjb. The wave equation can be solved as a series expansion 
in 1/q 2 . The zeroth order term is given by replacing F by F (eq. tH}). For A < |, F (z) has a smooth limit (albeit 
unphysical) as b — > oo (eq. After rescaling and letting 6 -» oo, we obtain the wave equation 



^4" + \z 



2A-2 



w 2 LA = 



(60) 



where we used F(z/b) — > F(0) = 1, as b — > oo. For 1 < A < |, there are two linearly independent solutions, A±, 
distinguished by their asymptotic behavior, 



4- 1 » A 



A± ~ e 

The general solution can be written as a linear combination, 

.4 = c + A+ + c~A 
Applying the boundary condition (|59")l . we deduce 



z ^ oo 



(61) 
(62) 
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so at very low temperatures, 

c+ = (64) 

i.e., A — > as z — > oo. 

For A = |, we obtain the exact explicit limit solution 

A{z) = A- (z) = Ai(bz - Co 2 ) (65) 

whereas A+(z) = Bi(&z — Co 2 ), with arbitrary normalization, where we restored the scaling parameter b. 

In this (unphysical) limit, the quasinormal frequencies have moved to the real axis yielding an infinite set of normal 
frequencies which are solutions of 

Ai(-w 2 ) = (66) 

Thus we obtain an infinite tower of real frequencies given by the zeroes of the Airy function. 
We deduce the limit conductivity (as q,b —> oo) 

i A^-Co 2 ) 

Ul Al(— Ul z ) 

The real frequencies that solve (j6"6")l are the poles of the conductivity. Notice that Sftcr = 0, except at the poles of Grx 
where -ftcr diverges as a (5-function. 

This expression is unphysical, as we have already explained (the probe limit breaks down in the limit b — > oo), but 
it is useful for computational purposes. We shall use it as a starting point to calculate the conductivity of a physical 
system at low temperatures. 

At low temperatures, we can calculate the first-order correction analytically by considering the b — > oo wave equation 
(1601) as the zeroth-order equation. Then for the first-order correction SA to the potential at low temperatures, we 
obtain from (j5"8")) . 

- SA" + \z - Co 2 ]6A = - — 1 .. HiA (68) 



3(1 -z 3 ) 



where 



d z 2 (l — 15z — 12z 2 — 10z 3 ) 

Ui = z \9z - 2(1 + z + z 2 )iu>] — + 3fe 3 z(2Fi(z) + z 3 ) - (1 + 2z + 3z 2 )iuj + — ^ -to 2 (69) 

dz 3(1 + z + z z ) 

The first-order potential leads to quasinormal modes which are zeroes of A + SA. Thus the real frequencies (p6l of 
the (unphysical) b — > oo limit get shifted at low temperatures away from the real axis. We obtain lj — > Co + 5Co, where 

7rBi(-w 2 ) f 1 dz A . /k _ a2vi/ 2 



8& = 3^/(^2) J Q Y^kiibz-Co^Mibz-Co 2 ) (70) 

This first-order expression is valid for low frequencies. As we heat up the system, most modes disappear and we are 
left with a finite number of quasinormal modes. Their number decreases as we increase the temperature. Conversely, 
as we cool down the system, (|70|) becomes increasingly accurate for an increasing number of modes. These modes 
shift toward the real axis (8b) — > 0) as we lower the temperature. In the low temperature regime ([2"T]) . the large q is, 
the more spikes one obtains for <r(w). 

This shifting of quasinormal modes can be seen in plots of the conductivity. As the mode frequency approaches the 
real axis, the corresponding spike in the plot of the imaginary part of the conductivity becomes more pronounced. To 
demonstrate this, we calculated the conductivity using the first-order approximation (|32[) to the scalar field. In figure 
we show the imaginary part of the conductivity at temperature T/T c w .1 and compare with the exact numerical 
solution. The agreement is very good even at such high temperature at which only one quasinormal mode is left. 
Unfortunately, this is the low temperature limit attained by numerical analysis as numerical instabilities prohibit one 
from lowering the temperature further. Using our analytic results, we see in figure [5] the emergence of an increasing 
number of poles as we lower the temperature to T/T c ss .06 and .04. Finally in figure [9] we compare the lower 
temperature (T/T c ~ .01) result with the b — > oo limit analytic expression (|67[) . thus demonstrating convergence. 

For A > 3/2, the potential is 

V = 6 2A z 2A ~ 2 (l - z 3 )F(bz) (71) 

with F given approximately by (|4"7| . It attains a maximum of order 6 2 ( 2 ~ A ) for A < 2. Therefore, at zero temper- 
ature it has infinite height. However, the width becomes infinitely narrow leading to a finite tower of poles for the 
conductivity (quasinormal modes). In the zero temperature limit, the number of modes increases as one approaches 
the BF bound and decreases away from it. For A > 2, the height of the potential becomes finite at zero temperature. 
It turns out that the potential is too narrow to possess bound states, so no poles exist for A > 2. 
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FIG. 7. The imaginary part of the conductivity in d = 3 using the expression (I32|l for the scalar field (dotted line) compared 
with the exact numerical solution (solid line) at ^- « .1 
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FIG. 8. The imaginary part of the conductivity vs. frequency in d = 3 using the expression (32) for F at 2- « .05 (left), .04 
(right). As the temperature decreases, poles move toward the real axis. 



B. Four dimensions 



The d — 4 case is similar. Working as in the d — 3 case, in the (unphysical) b — > oo limit, the wave equation for 
A = 2 (at the BF bound) in the probe limit reduces to 



A" - - A' - [b A z 2 - u 2 ]A = 



(72) 



whose acceptable solution can be written in terms of a Whittaker function, 

A = L (b 2 z 2 ) 



(73) 



(The other solution diverges asz-> oo.) At the boundary (z — > 0), it has a logarithmic divergence which we need to 
subtract before we can calculate QNMs and the conductivity |4|]. The conductivity is then given by 



2 A 2 iu 
a{oj) = — — + — 
iu A 2 



(74) 



where 



A{z) = A + A 2 b 2 z 2 - A ^b 2 z 2 \n(b 2 z 2 ) + 



(75) 
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FIG. 9. Comparison of the imaginary part of the conductivity in d = 3 using the expression (J32J) for F at ^ R! .01 (dotted 
line) and the b — > oo limit (|67[l (solid line) . 
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FIG. 10. The imaginary part of the conductivity in d — 4 using the expression (J32J for the scalar field (dotted line) compared 

T 



with the exact numerical solution (solid line) at ^- « .17 



with an arbitrarily chosen cutoff. 

Using the expansion for small arguments, 



i(&V) 



we deduce the limit conductivity as q, b — > oo, 



jl - [1 + 0? (2 7 - 1 + ln(6 2 z 2 ) + ^(1 - ^ 2 /4))] ^ j + 



1 



<t(w) = \- iuj 

iuj 



2 7 --+^(l-^ 2 /4) 



(76) 



(77) 



We have a pole at uj = 0, as expected and an infinite tower of real poles determined by the poles of the digamma 
function. The poles have real frequencies 



= 2V^ , n= 0,1,2,. 



(78) 



These poles are unphysical, but are good approximations to the quasinormal modes at low temperatures, in the 
physical regime (|21[) . As we increase the temperature, these modes move away from the real axis. At any given 
temperature we have a finite number of such modes. As we lower the temperature, their number increases and the 
quasinormal modes approach the poles (1781) on the real axis. To demonstrate this, we have calculated the conductivity 
using the first-order approximation (I32[) to the scalar field. In figure [TU] we compare with numerical results at 
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FIG. 11. The imaginary part of the conductivity vs. frequency in d = 4 using the expression (J32J) for F at « .1 (left), .04 
(right). 



temperature T/T c « .17 and find good agreement. As we go to lower temperature, numerical instabilities arise and it 
is no longer possible to compare our analytic results with their numerical counterparts. We find convergence to the 
b — > oo limit (|77|) . but much slower than in d = 3. In figure fill we show the imaginary part of the conductivity at 
T/T c rj 0.1 and 0.04. As we lower the temperature, the number of poles increases and the poles shift to the right 
on the real axis approaching the limit values (|78[) which correspond to the b — > oo limit of the conductivity shown in 
figure [H 




0.5 1.0 1.5 2.0 2.5 3.0 3.5 



V<o 2 > 

FIG. 12. The imaginary part of the conductivity in the (unphysical) limit b — > oo in d — 4 (|77[1 . 



VI. CONCLUSION 

We discussed holographic superconductors in the probe limit when the scalar hair of the dual black hole is near the 
BF bound. Using the analytic tools developed in [15j, we analyzed the low temperature regime trying to go beyond 
the range of applicability of numerical techniques which fail at very low temperatures. Thus undeterred by numerical 
instabilities, we found that at low temperatures the condensate diverges as | In T| 5 where S depends on the dimension 
of spacetime (eq. (|46[) ) . This signals the breakdown of the probe limit at low temperatures even at the BF bound. The 
divergence is very mild which explains why it was missed in earlier numerical analyses. Even though the probe limit 
at zero temperature cannot be a physical state for scaling dimensions A < s, it is still useful to analyze it, because 
it is a good approximation to physical states at low temperatures and thus a good starting point for computational 
purposes. 

We calculated the low temperature conductivity at the BF found in the probe limit and found exact limit analytic 
expressions in d = 3, 4. Thus, we showed that as we lower the temperature, the conductivity has an increasing number 
of spikes (quasinormal modes) which approach the infinite tower of real poles determined by the zeroes of the Airy 
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function in d — 3 (eq. (IBB"!) 1 ) and the Gamma function in d = 4 (eq. (|78p). As we heat up the system, the number of 
spikes decreases and their positions move off of the real axis. 

The probe limit we studied here can be used as a zeroth-order contribution to a perturbative expansion in 1/q 2 . 
Our results can be extended in a systematic way to include back reaction to the bulk metric in order to analyze 
physical states. This will greatly facilitate the probe of the zero temperature limit and shed light on the ground state 
of the system, as numerical methods fail due to numerical instabilities. Work in this direction is in progress. 
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